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We  present  and  analyze  a  method  of  laser  frequency  stabilization  via  steady-state  patterns  of 
spectral  holes  in  Eu^^:Y2Si05.  Three  regions  of  spectral  holes  are  created,  spaced  in  frequency  by 
the  ground  state  hyperfine  splittings  of  The  absorption  pattern  is  shown  not  to  degrade 

after  days  of  laser  frequency  stabilization.  An  optical  frequency  comparison  of  a  laser  locked  to  such 
a  steady-state  spectral-hole  pattern  with  an  independent  cavity-stabilized  laser  and  a  Yb  optical 
lattice  clock  demonstrates  a  spectral-hole  fractional  frequency  instability  of  1.0  x  10  T  that 
averages  to  g  x  10“^^  at  r  =  73  s.  Residual  amplitude  modulation  at  the  frequency  of  the  RE 

drive  applied  to  the  fiber-coupled  electro-optic  modulator  is  reduced  to  less  than  1  x  10 fractional 
amplitude  modulation  at  r  >  1  s  by  an  active  servo.  The  contribution  of  residual  amplitude 
modulation  to  the  laser  frequency  instability  is  further  reduced  by  digital  division  of  the  transmission 
and  incident  photodetector  signals  to  less  than  1  x  10“^®  at  r  >  1  s. 


Frequency  stable  laser  local  oscillators  (LLOs)  are  key 
tools  in  the  field  of  metrology.  Applications  of  such  LLOs 
include  optical  atomic  clocks  [1-3],  tests  of  general  rel¬ 
ativity  [4,  5],  searches  for  variation  of  fundamental  con¬ 
stants  [6,  7],  and  relativistic  geodesy  [8,  9].  The  best 
stabilized  lasers  to  date  [10,  11]  are  obtained  by  locking 
their  frequencies  to  Fabry-Perot  reference  cavities,  and 
their  stability  is  intrinscially  limited  by  thermomechani¬ 
cal  length  fluctuations  of  the  cavity  [12,  13]. 

Higher  stability  might  be  possible  by  the  method  of 
spectral- hole  burning  (SHE)  [14-16],  where  a  narrow  op¬ 
tical  transition  of  rare-earth  ions  doped  into  a  cryogeni- 
cally  cooled  crystal  is  used  as  the  frequency  reference. 
Strains  in  the  crystalline  host  inhomogenously  broaden 
the  narrow  absorption  line  by  several  gigahertz.  For  laser 
frequency  stabilization,  narrow  transparencies  (spectral 
holes)  are  created  within  the  absorption  line  by  frequency 
selective  optical  pumping  [17].  These  narrow  spectral 
holes  can  be  used  for  laser  frequency  stabilization  by  ac¬ 
tively  steering  the  frequency  of  the  laser  to  maximize  the 
transmission  through  the  crystal. 

Spectral  holes  in  Eu^+:Y2Si05  [18]  are  particularly 
promising  for  laser  frequency  stabilization.  At  4  K, 
this  material  supports  spectral  holes  at  580  nm  with 
linewidths  as  narrow  as  122  Hz  [19]  and  lifetimes  of  10^ 
s  [18].  The  frequency  shifts  due  to  fluctuations  in  the 
ambient  magnetic  and  electric  fields,  temperature,  pres¬ 
sure,  and  acceleration  are  all  small  enough  to  allow  laser 
frequency  instabilities  at  the  10“^^  fractional  frequency 
level  [20,  21].  However,  prior  laser  frequency  stabiliza¬ 
tion  experiments  with  Eu^+:Y2Si05  have  been  limited 
to  run-times  of  a  few  thousand  seconds  due  to  degrada¬ 
tion  of  the  spectral  holes  caused  by  noise  and  saturation 
of  the  probe  laser  [22]. 

In  this  work,  we  demonstrate  frequency  stabilization  to 
a  steady-state  pattern  of  spectral  holes  in  1.0  atomic  % 
Eu^+:Y2Si05  crystals  at  temperatures  near  4  K.  This 
pattern  consists  of  three  regions  of  spectral  holes  spaced 
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FIG.  1.  Absorption  of  a  steady- state  spectral  hole  pattern 
with  three  regions  of  spectral  holes.  The  left  inset  shows  a 
level  diagram  of  ^^^Eu^^  positioned  at  crystallographic  site  1. 
The  right  inset  shows  the  absorption  (blue)  and  corresponding 
PDH  error  signal  (red)  of  a  subset  of  the  center  region. 


in  frequency  by  46.2  MHz  and  34.6  MHz  (see  Figure  1), 
corresponding  to  the  ground-state  hyperfine  splittings  of 
^^^Eu^+  positioned  at  crystallographic  site  1  [23,  24]. 
The  Eu^+  population  reaches  steady-state  as  the  spectral 
holes  in  the  three  regions  are  burned.  Additional  inter¬ 
leaved  probing  does  not  modify  the  absorption  spectrum 
because,  for  a  subset  of  the  ^^^Eu^+  ions  with  certain  de¬ 
tunings  within  the  inhomogeneously  broadened  absorp¬ 
tion  line,  all  three  of  the  hyperfine  ground  states  are  ad¬ 
dressed  by  the  laser  as  it  probes  spectral  holes  in  the 
three  regions.  Absorption  does  not  accumulate  between 
the  spectral  holes  because  they  are  spaced  by  only  a  few 
times  the  minimum  observed  spectral-hole  width.  This 
persistent  and  self-regenerating  pattern  is  achieved  by 
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FIG.  2.  Schematic  diagram  of  the  spectral- hole  burning  ex¬ 
periment.  Solid  yellow  lines  denote  580  nm  laser  propaga¬ 
tion,  and  dashed  red  lines  denote  1160  nm  laser  propagation. 
AOM:  acousto-optic  modulator,  f-EOM:  fiber-coupled  electro¬ 
optic  modulator,  DET:  detector,  ECDL:  external  cavity  diode 
laser. 


EIG.  4.  Erequency  feedback  and  residual  amplitude  modu¬ 
lation  servo  schematic  diagram.  Black  arrows  denote  elec¬ 
tronic  signal  propagation,  and  the  dashed  box  denotes  cal¬ 
culations  performed  digitally  on  a  field-programmable  gate 
array  (EPGA).  LP:  low-pass  filter,  HP:  high-pass  filter,  LO: 
local  oscillator,  NE:  notch  filter. 


EIG.  3.  Pulse  sequence  of  the  interleaved  probe  scheme  for 
100%  duty  cycle  laser  frequency  stabilization.  Erequency  feed¬ 
back  from  the  previous  probe  pulse  on  each  crystal  is  applied 
to  the  pre-stabilization  cavity  AOM  at  the  times  indicated  by 
vertical  arrows. 


utilizing  probe  pulses  whose  center  frequency  coincides 
with  an  absorption  minimum.  Weaker  phase  modulation 
sidebands  are  spaced  to  coincide  with  absorption  maxima 
to  yield  a  Pound-Drever-Hall  type  error  signal.  Using 
this  spectral-hole  pattern,  laser  frequency  stabilization 
experiments  can  be  made  to  run  indefinitely. 

Our  SHB  LLO  employs  a  two-stage  laser  frequency 
stabilization  scheme,  as  shown  in  Figure  2  and  similar 
to  that  described  in  Reference  [22].  A  1160  nm  exter¬ 
nal  cavity  diode  laser  (ECDL)  is  frequency  doubled  to 
580  nm,  and  pre-stabilized  to  a  Fabry-Perot  cavity  [25]. 
The  pre-stabilized  laser  fractional  frequency  instability 
is  3  X  10“^^  for  0.1  s  <  r  <  1  s.  This  light  is  split  into 
two  beams  that  are  sent  through  two  independent  fre¬ 
quency  and  intensity  tuneable  acousto-optic  modulators 
(AOMs)  in  order  to  burn  and  probe  spectral  holes  in  two 
Eu^+:Y2Si05  crystals.  A  laser  intensity  of  2  juw/cm^ 
incident  on  the  crystals  is  used  for  both  burning  and 


probing  for  all  of  the  experiments  discussed  in  this  letter. 
Two  fiber-coupled  electro-optic  modulators  (f-EOMs)  ap¬ 
ply  sidebands  to  the  beams  at  15  kHz  for  crystal  one  and 
21  kHz  for  crystal  two,  for  generation  of  a  Pound-Drever- 
Hall  (PDH)  error  signal  [17,  26].  Detectors  (DETs)  II 
and  12  measure  the  power  incident  on  the  two  crystals, 
which  is  stabilized  by  feedback  to  the  RF  power  of  the 
probe  AOM  drives.  Detectors  T1  and  T2  measure  the 
probe  light  transmitted  through  the  crystals.  For  mea¬ 
surements  of  the  laser  frequency  stability,  visible  laser 
light  is  sent  to  an  independent  high-finesse  Fabry-Perot 
reference  cavity  [27],  and  infrared  laser  light  is  sent  to  a 
frequency  comb  for  comparison  against  other  frequency 
references,  each  through  a  Doppler-cancelled  optical  fiber 
[28]. 

During  a  SHB  laser  frequency  stabilization  run,  sev¬ 
eral  hundred  spectral  holes  are  created.  A  hole  is  then 
randomly  selected,  and  the  laser  spends  a  few  millisec¬ 
onds  probing  the  center  of  the  spectral  hole.  Frequency 
corrections  from  the  two  crystals  are  averaged  and  ap¬ 
plied  to  the  pre- stabilization  cavity  AOM.  In  order  to 
reduce  laser  frequency  degradation  due  to  the  Dick  ef¬ 
fect  [29] ,  probe  pulses  for  the  two  crystals  are  interleaved 
with  a  100%  duty  cycle,  as  shown  in  Figure  3.  To  reduce 
frequency  instability  caused  by  error  signal  transients  at 
the  beginning  of  each  probe  pulse,  the  error  signal  gen¬ 
erated  during  the  first  several  hundred  microseconds  of 
each  pulse  is  not  used  for  frequency  feedback. 

A  schematic  diagram  of  the  error  signal  generation  and 
residual  amplitude  modulation  (RAM)  servo  is  shown 
in  Figure  4.  The  transmission  detector  signal  is  digi¬ 
tally  divided  by  the  incident  detector  signal  on  a  field- 
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programmable  gate  array  (FPGA)  in  order  to  suppress 
laser  intensity  noise.  This  also  provides  some  suppression 
of  the  RAM  caused  by  the  f-EOM.  After  filtering,  this 
signal  is  multiplied  with  the  local  oscillator  to  produce 
the  PDH  error  signal.  The  incident  detector  signal  is  in¬ 
dependently  demodulated  in  order  to  control  the  RAM 
through  a  servo- loop  by  applying  a  DC  bias  to  the  f-EOM 
[30].  This  suppresses  the  contribution  of  RAM,  converted 
to  fractional  frequency  instability  of  a  laser  locked  to  a 
steady-state  hole  pattern,  to  less  than  1  x  10“^^  at  aver¬ 
aging  times  greater  than  1  s,  and  to  less  than  1  x  10“^^ 
at  averaging  times  greater  than  100  s. 

To  reduce  frequency  instability  due  to  temperature 
fluctuations  [21],  the  Eu^+:Y2Si05  crystals  are  housed 
in  a  sealed  copper  pillbox  with  a  fixed  amount  of  he¬ 
lium  gas.  The  temperature  of  the  cryocooler  is  servo 
controlled  to  keep  the  copper  pillbox  within  5  mK  of  the 
first-order  insensitive  temperature  (3.941  K).  The  mea¬ 
sured  out-of-loop  temperature  instability,  converted  to 
fractional  frequency  instability,  is  less  than  1  x  10“^^  at 
averaging  times  greater  than  1  s,  and  less  than  1  x  10“^^ 
at  averaging  times  between  100  s  and  10^  s. 

The  time-evolution  of  an  unperturbed  spectral  hole  at 
3.941  K  is  shown  in  Eigure  5.  The  top  plot  of  Eigure 
5  shows  a  decrease  in  hole  depth  (in  units  of  aL  from 
maximum  absorption  of  the  crystal)  and  an  increase  in 
hole  width.  The  bottom  plot  of  Eigure  5  shows  the  time- 
evolution  of  the  normalized  area  of  the  spectral-hole.  The 
hole  area  is  fit  to  a  three-part  exponential  decay  with 
empirically  determined  time  constants  of  2.7  h,  53.3  h, 
and  856.2  h.  During  a  steady- state  laser  frequency  sta¬ 
bilization  run,  the  entire  spectral-hole  pattern  must  be 
probed  frequently  enough  to  re-burn  all  of  the  spectral 
holes  before  they  decay.  Spectral  diffusion  [31,  32]  causes 
the  spectral  holes  to  broaden,  limiting  the  minimum  fre¬ 
quency  spacing  of  holes  within  the  steady-state  pattern. 

The  absorption  pattern  of  the  three-region  steady-state 
spectral-hole  pattern  after  three  hours  of  laser  frequency 
stabilization  run-time  is  shown  in  Eigure  1.  Each  of  the 
three  regions  has  201  holes  spaced  by  6  kHz.  As  the  pat¬ 
tern  reaches  steady-state,  the  spectral  holes  increase  from 
an  initial  1.7  kHz  full-width  at  half- maximum  (EWHM) 
to  their  steady-state  EWHM  of  3.3  kHz  after  ten  hours 
of  run-time. 

In  contrast,  a  subset  of  the  single-region  spectral- hole 
pattern  used  in  previous  work  [22]  is  shown  in  Eig¬ 
ure  6,  after  three  hours  of  run-time.  The  603  holes 
of  this  pattern  are  spaced  by  270  kHz.  Starting  from 
an  initial  width  similar  to  that  of  the  three-region  pat¬ 
tern,  the  spectral  holes  have  broadened  to  5.6  kHz  after 
three  hours  of  run-time.  The  spectral-holes  continue  to 
broaden  indefinitely,  which  degrades  the  frequency  stabil¬ 
ity  of  the  laser  locked  to  this  single-region  spectral-hole 
pattern. 

Interspersed  periodically  between  probe  pulses  used  for 
laser  frequency  stabilization  are  scan  pulses  of  individual 


FIG.  5.  Time  evolution  of  a  spectral  hole  in  the  absence  of 
degradation  due  to  probing.  The  top  plot  shows  a  decrease 
in  the  absorption  depth  (red  circles,  left  axis)  and  an  increase 
in  the  full-width  at  half-maximum  (FWHM)  (blue  triangles, 
right  axis)  with  a  fit  of  the  width  to  a  function  A  B(1  — 
g-t/r)  bottom  plot  shows  the  area  of  the  spectral  hole 

normalized  to  the  initial  area.  The  data  is  fit  to  a  three-part 
exponential  decay. 


FIG.  6.  Absorption  (blue)  and  corresponding  PDH  error  sig¬ 
nal  (red)  of  two  out  of  a  total  of  603  spectral  holes  in  a  single¬ 
region  spectral-hole  pattern. 

spectral  holes,  which  allow  for  the  hole  shape  and  error 
signal  slope  to  be  monitored  as  a  run  progresses.  The 
top  plot  of  Figure  7  shows  a  comparison  between  the 
slope  of  the  error  signal  of  a  single-region  pattern  and 
a  three-region  steady-state  pattern,  as  a  function  of  the 
duration  of  time  that  the  pattern  has  been  used  for  laser 
frequency  stabilization.  These  values  are  an  average  of 
the  error  signal  slope  at  the  center  of  each  spectral  hole 
within  the  pattern.  For  both  patterns,  the  error  signal 
slope  reaches  a  maximum  as  the  absorption  at  the  center 
of  the  holes  approaches  zero,  then  decreases  as  the  holes 
start  to  broaden.  The  three-region  steady-state  pattern 
asymptotically  approaches  a  non-zero  error  signal  slope, 
with  a  time  constant  that  indicates  how  long  it  takes  the 
pattern  to  reach  steady-state.  The  single-region  pattern 
error  signal  slope  decreases  until  the  lock  can  no  longer 
be  maintained  after  approximately  70  hours  of  run-time. 
The  bottom  plot  of  Figure  7  shows  the  1  second  averaging 
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FIG.  7.  Comparison  of  a  single-region  pattern  (blue  squares) 
versus  a  three-region  steady-state  spectral-hole  pattern  (red 
circles).  The  top  plot  shows  the  error  signal  slope  averaged 
over  all  of  the  holes  in  the  pattern,  and  the  bottom  is  the 
1  second  crystal  versus  crystal  Allan  deviation  with  linear 
frequency  drift  removed.  Both  are  plotted  with  respect  to 
the  run-time  of  the  laser  frequency  stabilization. 


time  Allan  deviation  of  the  differential  frequency  noise 
of  the  two  laser  beamlines  independently  locked  to  the 
two  crystals,  after  linear  drift  has  been  removed.  The 
Allan  deviation  values  have  been  divided  by  v/2,  which 
is  an  upper  bound  on  the  cystal  with  better  spectral-hole 
frequency  stability.  After  the  three-region  pattern  has 
reached  steady-state,  the  1  second  differential  fractional 
frequency  instability  never  rises  above  2  x  10“^^,  out  to 
a  run-time  of  almost  ten  days.  The  differential  frequency 
instability  of  the  single-region  pattern  steadily  increases, 
rising  above  1  x  10“^^  at  33  hours  of  run-time. 

We  have  measured  the  frequency  stability  of  a  laser 
locked  to  a  three-region  steady-state  pattern  of  spectral 
holes  in  Eu^+:Y2Si05  using  a  three-cornered  hat  method 
[33] .  The  stability  is  measured  relative  to  an  independent 
Fabry-Perot  reference  cavity  [27]  and  a  Yb  optical  lattice 
clock  with  a  stability  of  3.2  x  10“^^  [34],  xhe  frac¬ 

tional  frequency  instability  of  the  laser  stabilized  to  the 
spectral  holes  is  1.0  x  10“^^  for  0.01  s  <  r  <  20  s, 

and  reaches  x  10“^^  at  r  =  73  s,  as  shown  in  Fig¬ 

ure  8.  The  duration  of  this  measurement  was  limited  by 
the  duration  of  time  the  frequency  comb  remained  mode- 
locked.  The  crystal  versus  crystal  fractional  frequency 
instability  averages  to  4. big  7  x  10“^^  at  r  =  2000  s  with 
a  0.5  mHz/s  linear  drift  removed.  This  Allan  deviation 
of  the  frequency  difference  between  the  two  crystals  has 
been  divided  by  \/2  to  put  an  upper  bound  on  the  fre¬ 
quency  stability  of  the  better  of  the  two  crystals.  Figure 
8  also  shows  the  estimated  contribution  of  temperature 
instability,  Doppler  shifts,  RAM,  and  detection  noise  to 
the  frequency  instability  of  the  SHB  LLO.  Temperature 
instability  is  measured  using  an  out-of-loop  sensor  and 
converted  to  frequency  instability  assuming  operation  at 
5  mK  away  from  the  first-order  insensitive  temperature. 


FIG.  8.  Fractional  frequency  instability  of  a  laser  locked  to  a 
steady-state  pattern  of  spectral  holes  in  Eu^^:Y2Si05  .  Ab¬ 
solute  laser  frequency  instability  from  a  three-cornered  hat 
measurement  is  shown  (blue  squares),  as  well  as  the  differen¬ 
tial  instability  of  a  crystal  versus  crystal  frequency  compari¬ 
son  (triangles)  where  much  of  the  technical  noise  is  common¬ 
mode.  The  instability  of  the  crystal  versus  crystal  measure¬ 
ment  is  worse  at  small  averaging  times  because  the  duty-cycle 
is  50%,  while  the  absolute  spectral  hole  stability  measurement 
uses  feedback  from  both  crystals  to  achieve  100%  duty  cycle. 
The  expected  frequency  noise  due  to  temperature  fluctuations 
(dashed  grey),  Doppler  shifts  due  to  motion  of  the  pillbox 
(solid  grey),  as  well  as  RAM  and  detection  noise  (open  red 
circles)  are  also  shown. 


Doppler  shifts  introduced  by  the  motion  of  the  cryostat 
are  measured  interferometrically  by  retroreflecting  a  laser 
off  of  a  mirror  mounted  on  the  copper  pillbox.  These 
sources  of  laser  frequency  instability  are  shown  not  to 
limit  the  present  performance  of  this  spectral-hole  burn¬ 
ing  experiment. 

We  have  demonstrated  steady-state  SHB  laser  fre¬ 
quency  stabilization  with  a  short-term  fractional  fre¬ 
quency  instability  ofl.0xl0“^^r“^/^  that  averages  down 
to  8.5li;g  X  10“^^  at  r  =  73  s.  The  steady-state  spectral- 
hole  pattern  successfully  prevents  degradation  of  the  laser 
frequency  stability  out  to  run-times  of  many  days.  RAM 
is  suppressed  to  less  than  1  x  10“^  fractionally  at  1  s  aver¬ 
aging  time  using  an  active  servo,  and  the  digital  division 
of  transmitted  and  incident  photodetector  signals  further 
suppresses  laser  frequency  instability  due  to  RAM. 
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